Dynamic user equilibrium with elastic demand and bounded 
rationality: Formulation, qualitative analysis and computation 
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Abstract 

This paper is concerned with two extensions of the simultaneous route- and- departure 



choice (SRDC) dynamic user equilibrium (DUE) originally proposed by Friesz et al. ( 1993 1. 
The first extension is an SRDC DUE with elastic travel demand, denoted by E-DUE, in 
which the trip matrix is determined endogenously. The second extension is a bounded 
rationality SRDC DUE, denoted by BR-DUE, in which travelers do not always seek the 
least costly route-and-departure-time choice. We propose to analyze and compute these 
two types of equilibrium models using the emerging mathematical paradigm of differen- 



tial variational inequality (Pang and Stewart 20081. As we shall explain, both extensions 



can be expressed in continuous time as a variational inequality (VI) and a differential 
variational inequality (DVI). The VI is a canonical mathematical form which allows ex- 



istence results to be established by means of Brouwer's fixed point theorem (H an et al 
2013c|d| ). The DVI formulation enables us to reformulate both E-DUE and BR-DUE as 



fixed-point problems in proper Hilbert spaces, by employing necessary optimality condi- 
tions for optimal control problems. We devise two fixed-point iterative algorithms for 
computing E-DUE and BR-DUE in continuous time followed by some qualitative analysis 
on the convergence issue. Both algorithms can accommodate a wide range of analytical or 
simulation-based dynamic network loading submodels. We demonstrate the efficiency and 
convergence of these algorithms with several numerical studies on the Sioux Falls network. 
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1 Introductory Remarks 

This paper is concerned with two extensions of the simultaneous route-and-departure choice 



dynamic user equilibrium (SRDC DUE) originally proposed in Friesz et al. (1993) and dis- 



cussed subsequently by Friesz and Mookherjee (2006); Friesz et al. (2011 2013); and Friesz 



and Meimand (2013). The first extension relaxes the assumption of fixed trip volumes by 



considering elastic travel demands among origin-destination pairs. As such, it is concerned 
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with SRDC-DUE for which travel cost, including delay as well as early/late arrival penalties, 
are equilibrated and demand is determined endogenously. The second extension incorporates 



the concept of bounded rationality (BR) proposed by Simon (1957, 1990 1991) in the mod 



eling of travel behavior. Bounded rationality-based traffic equilibrium models operate under 
the assumption that travelers, viewed as Nash agents, do not behave in a completely rational 
manner. 



1.1 Dynamic user equilibrium with elastic demand 

Most of the studies of DUE reported in the dynamic traffic assignment (DTA) literature are 
about dynamic user equilibrium with constant travel demand for each origin-destination pair. 
It is, of course, not generally true that travel demand is fixed, even for short time horizons. 



Arnott et al. (1993) and Yang and Huang (1997) directly consider elastic travel demand 



in the context of a single bottleneck. Yang and Meng (1998) extend a simple bottleneck 



model to a general queuing network with known elastic demand functions for each origin- 
destination (OD) pair. They employ a space-time expanded network (STEN) representation 
of the network loading submodel. Wie et al. (2002) study a version of the dynamic user 



equilibrium with elastic demand, using a complementarity formulation that requires path 
delays to be expressible in closed form. Szeto and Lo (2004) study dynamic user equilibrium 
with elastic travel demand when network loading is based on the cell transmission model 
(CTM); their formulation is discrete-time in nature and is expressed as a finite-dimensional 
variational inequality (VI). The VI is solved with a descent method under the assumption 



that the delay operator is co-coercive, or to have the Dunn property (Farouq, 2001). Han 



et al. (2011) study dynamic user equilibrium with elastic travel demand for a network with 



a single origin-destination pair whose traffic flow dynamics are also described by CTM; the 
CTM is chosen to accommodate the discrete-time complementarity formulation of the user 
equilibrium model. 



Although Friesz et al. (2011 ) show that analysis and computation of dynamic user equilib- 



rium with constant travel demand is tremendously simplified by stating it as a a differential 
variational inequality (DVI), they do not discuss how elastic demand may be accommodated 



within a DVI framework. Friesz and Meimand (2013) later extend the DVI formulation to an 



elastic demand setting, although that paper does not discuss computation of E-DUE, which is 
our main focus herein. Such a DVI formulation of elastic demand DUE is not straightforward. 
In particular, the DVI presented therein has both infinite-dimensional and finite-dimensional 
terms. Moreover, for any given origin-destination pair, inverse travel demand corresponding 
to a dynamic user equilibrium depends on the terminal value of a state variable representing 
cumulative departures. The DVI formulation achieved in that paper is significant because it 
allows the still emerging theory of differential variational inequalities to be employed for the 
analysis and computation of solutions of the elastic-demand DUE problem when simultaneous 
departure time and route choice are within the purview of users, all of which constitutes a 
foundation problem within the field of dynamic traffic assignment. Han et al. (2013d) study 



the existence of E-DUE proposed in Friesz and Meimand (2013) by formulating it as an 



infinite-dimensional variational inequality in a Hilbert space, and by applying an extension of 



Brouwer's fixed theorem (Browder, 1968). 



A good review of recent insights into abstract differential variational inequality theory, 
including computational methods for solving such problems, is provided by |Pang and Stewart| 
(2008). Also, differential variational inequalities involving the kind of explicit, agent-specific 



control variables employed herein are presented in Friesz (2010). 



One of the main goals of this paper is to establish the most general framework for com- 
puting simultaneous route-and-departure choice E-DUE, regardless of the underlying traffic 
flow model or the dynamic network loading (DNL) subproblem employed. We also wish to 
provide analysis on the convergence of the proposed algorithm. To this end, we present and 
prove an equivalent fixed-point-problem (FPP) formulation of the SRDC E-DUE, and devise 
a fixed-point iteration scheme for computing it in continuous time. As stated in the literature 
review above, most computational schemes rely on a specific type of DNL model. To the best 
knowledge of the authors, the proposed fixed-point algorithm is the first general computational 
method that is applicable with any DNL submodel and elastic demand function under mild 
assumptions. 



1.2 Bounded rationality dynamic user equilibrium 

In dynamic traffic assignment (DTA) problems, the modeling of travelers' route and departure 



time choices have been greatly influenced by Wardrop's first principle (Wardrop, 1952), which 



states that road users behave in a completely rational way and seek to minimize their own 

disutility by adjusting route and/or departure time choices. There are multiple means of 

expressing the dynamic notion of Wardropian user equilibrium, such as variational inequality 

Friesz et al. 1993; Han et al. 2013c|d Szeto and Lo 2004), differential variational inequality 



Friesz et al., 2001, 2011, 2013), and nonlinear complementarity problem (Han et al. 2011 



Wie et al. , 2002). While admitting a number of canonical mathematical representations, the 



notion of complete rationality user equilibrium may enjoy less support from the point of view 
of human behavior and from empirical evidence. That is, travelers may not always choose 
the departure time and route that yield the minimum travel cost. Such a situation could 
be caused by (1) imperfect travel information; and (2) certain "inertia" in decision-making. 
Moreover, empirical studies suggest that in reality, drivers do not always follow the least costly 



route-and-departure-time choice (Avineri and Prashker, 2004). 



As a relaxation of the otherwise restrictive Wardropian assumption, the notion of bounded 



rationality is proposed by Simon ( 1957 , 1990 , 1991 ) and introduced to traffic modeling by Mah 



massani and Chang (1987). In prose, the notion of bounded rationality postulates a range 



of acceptable travel costs that, when achieved, do not incentivize travelers to switch their 



departure time or route choice. Such a range is phrased by Mahmassani and Chang (1987) as 



"indifference band". The width of such band, usually denoted by e, is either derived through 



a behavioral study of road users or calibrated by empirical observation (Hu and Mahmas- 



sani . 



1997 Mahmassani and Jayakrishnan 1991 Mahmassani and Liu, 1999). In general 



e could depend on origin-destination pairs and/or travel commodities. Bounded rationality 
user equilibrium plays an important role in transportation modeling and particularly in static 



traffic assignment (STA), with an incomplete list of papers including Di et al. (2013); Gifford 



(2010) and Marsden et al. (2012). It is also investigated via simulations-based approaches 



and Checherita (2007); Han and Timmermans (2006); Khisty and Arslan (2005); Luo et al. 



Mahmassani and Liu, 1999, Mahmassani et al., 2005) 



Although the concept of bounded rationality have been used extensively in STA and 
simulation-based models, there have been very few analytical results regarding bounded ratio- 



nality in dynamic traffic assignment (DTA). Szeto (2003) and Szeto and Lo (2006) are the first 



to propose bounded rationality dynamic user equilibrium (BR-DUE), although it is phrased 



differently as tolerance-based dynamic user optimal by the authors. In Szeto and Lo (2006), a 



route-choice (RC) BR-DUE is formulated as a discrete-time nonlinear complementarity prob- 



lem. They propose a heuristic route-swapping algorithm adapted from Huang and Lam ( 2002 ) 



to solve the RC BR-DUE. Since then, no further extension of such BR-DUE has been made 
in the literature. In this paper, we propose a (differential) variational inequality framework 
in which the simultaneous-route-and-departure choice (SRDC) BR-DUE can be analyzed and 
computed; thus advancing the model of BR-DUE to accommodate more general travel choices. 
Moreover, our analysis and computational tool do not require specific forms of the dynamic 
network loading (DNL) submodel, thus can be applied to traffic modeling with a wide range 
of analytical or simulation-based DNL procedures. As minor extension of the proposed model, 
we also consider the cases where the "tolerance", or e, can depend on origin-destination pairs 
and travel commodities. All these extensions, as we shall illustrate in this paper, can be easily 
handled by the mathematical paradigms of variational inequality and differential variational 
inequality. 

1.3 The VI and DVI framework for analyzing DUE 

In this paper, we present a unified theory and a general framework for analyzing the aforemen- 



tioned two extensions of SRDC DUE proposed by Friesz et al.| { 1993 >, and many others to be 
yet identified. In particular, we present variational inequality (VI) and differential variational 
inequality (DVI) formulations of dynamic user equilibrium with elastic demand (E-DUE) and 
bounded rationality dynamic user equilibrium (BR-DUE), when simultaneous departure time 
and route choice are within the purview of users. As we shall demonstrate in the paper, 
the notions of elastic demand and bounded rationality can be easily captured by variational 
inequalities defined in appropriate Hilbert spaces. These VI formulations also facilitate the 
derivation of existence results in continuous time, by invoking a version of Brouwer's fixed 



point theorem (Browder, 1968), see also Han et al. (2013c) for an example. 



The computation of a DUE is best facilitated by the mathematical paradigm of differential 



variational inequality (DVI) (Pang and Stewart, 2008) and emerging computational algorithms 



associated therein. The VI formulation, either for E-DUE or BR-DUE, easily transfers to a 
DVI formulation by means of a two-point boundary value problem. This paper employs the 
optimal control framework to further establish fixed-point-problem (FPP) formulations of 
both E-DUE and BR-DUE. As far as we know, the proposed FPP formulations are wholly 
original in the literature and are established for almost all types of dynamic network loading 
submodels. As immediate consequences of the FPP formulations, two fixed-point iteration 
schemes are proposed that computes E-DUE and BR-DUE respectively in continuous time. 
Each intermediate iteration of the proposed algorithms is equivalent to a linear-quadratic 
optimal control problem, which we manage to solve in closed-form using necessary optimality 
conditions for optimal control problems. Qualitative analysis regarding existence of solutions 
as well as convergence of the iterative algorithms are provided as well. The proposed algorithms 
are tested numerically; the results demonstrate their effectiveness and solution quality. 
The main findings made in this paper include: 

• the expression of the simultaneous route-and-departure choice (SRDC) E-DUE problem 
as a fixed-point problem that subsumes several existing models we reviewed above; 

• a continuous-time algorithm for the above fixed-point problem that, unlike many exist- 
ing algorithms, accommodates a wide range of analytical or simulation-based dynamic 
network loading procedures; 

• the variational inequality and differential variational inequality representations of SRDC 
BR-DUE, which simplify the analysis of BR equilibrium such as existence, and allows 
heterogeneity in the indifference band to be captured; 



• the expression of the SRDC BR-DUE as a fixed-point problem and a continuous-time 
algorithm that computes SRDC BR-DUE, which also encapsulates almost all types of 
dynamic network loading models; and 

• convergence results for the two continuous-time algorithms mentioned above for a limited 
class of path delay operators. 

1.4 Organization 

The rest of this paper is organized as follows. Section [2] provides a general review of the nota- 
tions and assumptions employed in the notion of dynamic user equilibrium proposed by |Friesz 



et al.| (1993). Section [3] conducts detailed analysis of the simultaneous route-and-departure 



choice (SRDC) DUE with elastic demand (E-DUE) on a general network, using the theories 
of differential variational inequality and optimal control. A fixed-point formulation and al- 
gorithm are proposed that allow the E-DUE to be computed in continuous-time. In Section 
Hi we introduce, for the first time, the variational inequality (VI) formulation of SRDC DUE 
with bounded rationality (BR). As we explain, the notion of BR can be naturally captured by 
the VI formulation with appropriately defined principal operator. As an immediately conse- 
quence of the VI formulation, we present the differential variational inequality (DVI) and fixed 
point problem (FPP) formulations of BR-DUE. A fixed-point algorithm is also devised that 
computes BR-DUE in continuous time. Section [5] conduct qualitative analysis of both models 
by investigating the existence of solutions. Section [6] is concerned with the convergence of the 
proposed fixed-point algorithms. In Section [7j we conduct numerical tests on the fixed-point 
algorithms and present solution features as well as convergence results. 

2 Notation and Essential Background 

For the rest of this paper, the time interval of analysis is a single commuting period or "day" 
expressed as 

[to, t f ] c K 

where tf > to, and both to and tf are fixed. Here, as in all DUE modeling, the single most 
crucial ingredient is the path delay operator, which provides the delay on any path p per unit 
of flow departing from the origin of that path; it is denoted by 

D p (t, h) V P £V 

where V is the set of all paths employed by travelers, t denotes departure time, and h is a 
vector of departure rates. From these, we construct effective unit path delay operators ^f p (t, h) 
by adding the so-called schedule delay f(t + D p (t, h) — Ta); that is 

%(t, h) = D p (t, h) + f(t + D p (t,h)-T A ) V P £V (2.1) 

where Ta is the desired arrival time and Ta < tf. The function /(•) assesses a penalty 
whenever 

t + D p (t, h) £ T A (2.2) 

since t + D p (t, h) is the clock time at which departing traffic arrives at the destination of path 
p £ V. Note that, for convenience, Ta is assumed to be independent of destination. However, 
that assumption is easy to relax, and the consequent generalization of our model is a trivial 



extension. We interpret ^ p (t, h) as the perceived travel cost of drivers starting at time t on 
path p under travel conditions h. We stipulate that each 

%(-, h) ■ [t , t f ] — > K++ Vp g V 

is measurable, strictly positive, and square integrable. We employ the obvious notation 

*(•,/») = (%(-,h): P eV)e^l 

to express the complete vector of effective delay operators. 

It is now well known that path delay operators may be obtained from an embedded delay 
model, data combined with response surface methodology, or data combined with inverse 
modeling. Unfortunately, regardless of how derived, realistic path delay operators do not 
possess the desirable property of monotonicity; they may also be non-differentiable. 



2.1 Dynamic User Equilibrium Introduced by Friesz et al. (1993) 



For the comprehensiveness of our presentation, we recap in this section the notion of dynamic 
user equilibrium originally articulated by Friesz et al.| (|1993[) . 



We introduce the fixed trip matrix {Qij : (i, j) G W), where each Qij G K + is the 
fixed travel demand between origin-destination pair (i,j) G W, where W is the set of origin- 
destination pairs. Note that Qij represents traffic volume, not flow. Finally we let Vij CPbe 
the set of paths connecting origin-destination pair (i,j) G W. As mentioned earlier h is the 
vector of path flows h = (h p : p G V) ■ We denote the space of square integrable functions on 
the real interval [to, tf] by L 2 [to, tf]. We stipulate that each path flow is square integrable, 
that is 

h€(Ll[t ,t f ]f l 

where (L+[£o, tf]) is the positive cone of the |"P|-fold product of the Hilbert spaces L 2 [to, tf}. 
Each element of such set, h = (h p : p G V), is interpreted as a vector of departure-time 
densities, or more simply path flows, measured at the entrance of the first arc of the relevant 
path. It will be seen that these departure time densities are defined only up to a set of measure 
zero. With this in mind, let v denote a Lebesgue measure on [to, tf], and let V„(t) G [to, tf] 
denote the phrase for v-almost every t G [to, tf]. 
We write the flow conservation constraints as 

^2 / hp(t) dt = Q^ V(i, j) G W (2.3) 



P&V l: 



where (2.3 ) consists of Lebesgue integrals. Using the notation and concepts we have mentioned, 



the feasible region for DUE when the effective delay operators are known is 



An 



| h > : J^ I f h p (t) dt = Q^ V (i, j) G W \ C {L% [t , tf]) |P| (2.4) 

( pep* Jt0 ) 



In order to define an appropriate concept of minimum travel costs in the present context, we 
require the measure-theoretic analog of the infimum of a set of numbers. In particular, for 
any measurable function g : [to, tf] — > K, the essential infimum of g(-) on [to, tf] is given by 

essinf{#(,s) : s e [to, tf]} = sup{xG^: v{s G [to, tf] : g(s) < x} = 0} (2.5) 



Note that for each x > essinf{g(s) : s G [to, tf]} it must be true by definition that 

v{s G [t , t f ] : f{s) < x} > 

Let us define the essential infimum of effective travel delays, which depend on the path flows 
h: 

v p (h) = essinf {9 p (t, h) : t G [to, */]} > \/p G V (2.6) 

Vij {h) = mm {v p {h): peVij} V(i,j)eW (2.7) 



The following definition of dynamic user equilibrium is first articulated by Priesz et al 



(1993): 



Definition 2.1. (Dynamic user equilibrium) A vector of departure rates (path flows) 
h* G Ao is a dynamic user equilibrium if 

h;(t)>0, peVij =>- V p (t, h*) = Vij (h*) V„(t) g [to, t f ] (2.8) 

W^e denote this equilibrium by DUE(^>, Ao, [to,*/])- 



Using measure-theoretic arguments, Friesz et al. (1993) establish that a dynamic user 



equilibrium is equivalent to the following variational inequality under suitable regularity con- 
ditions: 

find h* G Ao such that 



Y^ / ^ P (t,h*)(h p -h* p )dt>0 

„,— n Jtn 



V/iG A 



► VI(9, A , [t , t/]) (2.9) 



3 Dynamic User Equilibrium with Elastic Demand (E-DUE) 

The general setup of DUE with elastic demand is similar to that of the fixed demand case, 
with the exception that the total travel demand Qij between an origin-destination (OD) pair 
(i, j) G W is no longer a prescribed constant. Rather, transportation demand is assumed to 
be expressed as the following invertible function 

Qij = Ftj[v] 

for each origin-destination pair (i, j) G W, where W is the set of all origin-destination pairs and 
v is a concatenation of origin-destination minimum travel cost Vij associated with (i, j) G W. 
That is, we have that 

v^ G K++, v = (vij : (i, j) G W) G ft™ 

Note that to say v^ is a minimum travel cost means it is the minimum cost for all departure 
time choices and all route choices pertinent to origin-destination pair (i,j) G W. Further note 
that Qij is the unknown cumulative travel demand between (i, j) G W that must ultimately 
arrive by time tf. 

We will also find it convenient to form the complete vector of travel demands by concate- 
nating the OD-specific travel demands to obtain 



Q = (Qij ■ (i, J) G W) G 3*+ 



|W| 



The inverse demand function for every (i,j) G W is 

Vij = @ij [Q] 

and we naturally define 

e = (Qij : (i, j) GW)efl 

As a consequence, we employ the following feasible set of departure flows when the travel 
demand between each origin-destination pair is unknown. 



A 



(h, Q) : h> 0, J^ / ftp(*)eft 

\V\ 



Qi 5 V(i,j)eW 



(3.10) 

where (L 2 [to, tf]) x 3?^+ is the direct product of the |"P|-fold product of Hilbert spaces 
consisting of square- integrable path flows, and the |W (-dimensional Euclidean space consisting 
of vectors of elastic travel demands. 

With preceding preparation, we are in a place where the simultaneous route-and-departure- 
choice dynamic user equilibrium with elastic demand can be rigorously defined, as follows. 

Definition 3.1. (Dynamic user equilibrium with elastic demand) A pair (h* , Q*) £ A, 
where h* is a vector of departure rates (path flows) and Q* is the associated vector of travel 
demands, is said to be a dynamic user equilibrium with elastic demand if for all (i, j) G W, 



h* p (t) > o, per, 



'.I 



* p (t,h*) = QijlQ*] 
V P (t,h*) > Q^Q*] 



V u (t) G [t , t f ] (3.11) 

V„(t) € [t , tf], VpePij (3.12) 



3.1 The Differential Variational Inequality Formulation of E-DUE 

This section is a close recap of the differential variational inequality (DVI) formulation of 



E-DUE problem proposed by Friesz and Meimand (2013). The notion of dynamic user equi- 



librium can be alternatively illustrated and analyzed using the mathematical paradigm of 



DVI (Pang and Stewart, 2008). In particular, as noted by Friesz et al. (2011), the demand 



satisfaction constraints in the presence of fixed demand can be easily rewritten as a two-point 
boundary value problem 



dy 



'j 



ell 



P&V i:j 



Vij(t ) = 0, yij(t f ) = Qij V(i, j) G W 



where yij (t) is interpreted as the cumulative departure associated with each origin-destination 
pair (i, j) G W by time t G [to, tf]. In the elastic demand case, we assume there are unknown 
terminal states yij(tf) for all (i,j) G W, which are the realized DUE travel demands. Moreover, 
for each origin-destination pair (i,j) G W, inverse travel demand is expressed as 



Vij = @ij[y(t f )] 



(3.13) 



where 



v = (wi(-): (i,j)ew) e(n%,t f ]) 



\w\ 



'H 1 ([to, tf]) is the Sobolev space consisting of weakly differentiable functions whose weak 
derivatives are square-integrable. Thus, the Vij(tf), for all (i,j) £ W, will be determined 



endogenously to the differential variational inequality presented subsequently in Theorem 3.2 



. Such an approach contrasts to the approach employed by Friesz et al. (2011 ) to study fixed 



demand DUE by making each yij{tf) an a priori fixed constant Qij. Accordingly, we introduce 
the following dynamics: 



dyij 
dt 






/-" 



yij{t ) = V(i,j)eW 



(3.14) 



Notice that (3.14) is an initial value problem without terminal constraints. As a consequence, 



we employ the following alternative form of the feasible set: 



Ai 



h>0: 



dyij 
dt 



Av\ 



J2 M*)> Vijfo) = V(i,i) £W)c(4 M/]) 1 ' (3.15) 

peVij 



Note that the feasible set Ai in (3.15) is expressed as a set of path flows since knowledge of h 



completely determines the state variables $/$,-, (i, j) £ W, that satisfy the inital value problem 
(3~14). 



Theorem 3.2. (E-DUE equivalent to a differential variational inequality) Assume 
ty p (-,h) : [to, tf] —> K++ is measurable and strictly positive for all p £ V and all h £ Ai. 
Also assume that the elastic travel demand function is invertible, with inverse @ij[Q] for all 
(i,j) £ W. A vector of departure rates (path flows) h* £ Ai is a dynamic user equilibrium 
with associated demand y*j(tf) if and only if h* solves the following DVI\9, 0, [to, */])■ 



find h* £ Ai such that V/i £ Ai, 

J2J% p (t,h*)(h p -h* p )dt- £ e y [y*(t/)] {yij{tf)-vUtf)) > o ^DVi{%e,[to,t t ]) 

(3.16) 

D 



p&v 



(ij)ew 



Proof. See Friesz and Meimand (2013). 



3.2 The fixed-point-problem reformulation of E-DUE 

In this section, we present and prove a fixed-point-problem (FPP) reformulation of E-DUE. 
Such a result is new in the literature and allows an iterative computational scheme to be 
devised based on the FPP. 



Since the differential variational inequality problem (3.16) has both infinite-dimensional 



and finite-dimensional terms, it is convenient for us to analyze this problem within an extended 
space that contains both (L 2 [to, tf]) and K' w ' as embedded subspaces. Therefore, we are 

prompted to define the product space E = (-L 2 [to, t/]) X W w ' which is a Hilbert space with 
the inherited inner product defined as 



1^1 f t f IW| 

{X, Y) E = Y, / &(*) • m{t) dt + Y, u i v i 

1=1 Jt ° j=l 

X = (60: • • ■ > C|P[(0i Wl) ■ • • i U |W|) S .B 

y = (r/i(-), • • • , V\v\(-), vi, ..., v m ) £ E 



(3.17) 



We also define the norm on the space E: 



1*11* = <*, x) 1 / 2 



' E 



Note that the elastic demand satisfaction conveyed by (3.10) can be trivially re- written as 
a two-point boundary problem, leading to the following alternative definition of the feasible 
set in the extended space E: 

Ao = ((/», Q) € (Ll[t , t f \f\ X »f I : ^ = £ h p (t), 

Wi(*o) = 0, yij{tf) = Q lj V(i,j)ew\ C E (3.18) 

Recalling the inverse demand function = (0y : (i, j) G W), we introduce the notation 

8- = (-9,, : (i, j) G W) : »f — ► ^L W| 
We also define the mapping 

T:Ao^E, (h,Q) >-> (V(-,h),e-(Q)) (3.19) 

where (h, Q) G A , *(•, /i) G (L+[*o, */])' W '' and @ _ (Q) G ^— • Sucn a mapping is clearly 
well-defined. With the preceding discussion, we are now ready to present and prove the main 
theorem in this section. 

Theorem 3.3. (E-DUE equivalent to a fixed-point problem) The fixed point problem 

X* = P~ ko [X* - aF{X*)] (3.20) 

is equivalent to the E-DUE problem, where X* = (h*, Q*) G {L 2 + [t , t/]) |p| X 3?!^', and P^[-] 
is the minimum norm projection onto Aq . a G K++ is a /ixed constant. 



Proof. The fixed-point problem (3.20) requires that 



X* = argmm{l\\X* - aT^X*) - X\\ 2 R : X £ A } (3.21) 



The right hand side of (3.21) amounts to the following optimal control problem when X* is 
given: 

min J(X) = r\Y, E [K(t)-a%(t,h*)-h p (t)] 2 dt 
X=[h,y{tf)) J to (i,j)eWpeVij 

+ \ E te(*/) -«e^ [/(</)] -Wi(*/)) 2 ( 3 - 22 ) 



2 
(<J)eW 



subject to 



EM*) V(t,i)€W (3.23) 



dt 



Vijfo) = V(i,j')€W (3.24) 

/i > (3.25) 
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The Hamiltonian for the above optimal control problem is 

H = 2 E E K(*)-«*p(*,^)-m*)] 2 + E A ^ E M*) 

(i,j)eW pev i:i (i,i)ew pePy 

for which the adjoint equations are 



r/A 



'./ 



<9# 



V(», j) 6W, p e T 7 ^, t e [to, */] 



(it dyij 

and the transversality conditions read: for all (i, j) G W, 

By the minimum principle, we enforce the following minimization problem 

mini? s.t. — h < 

for which the Kuhn- Tucker conditions are (we denote the K-T point by h**) 

- [h*(t) - a%(t, h*) - h**(t)] + Xij = Pp V(t, j) eW, pe p,y 

M;* = V(*,i)eW, peVi, 

P P > V(J,i)eW, pe^j 



(3.26) 



4(t / )+ae ij .[y*(t / )]+y ii (t / ) 

(3.27) 



(3.28) 
(3.29) 
(3.30) 



By virtue of ( |3.20| >, we have /i* = fe**. Eq uation ( |3.26| ) and ( |3.27[ ) then imply Ay 
aO~- [y*(t/)] • Therefore, we may re-state (3.28) as 



a%(t, h*) + a@r. [y^tf)] = p p V(«, j) € W, p£ 7>y 
Consequently, recalling 0~ = —0, we have by (3.29) that 

h* p (t)>0 =► p p = => * p (t,fc*) = By [^(t/)] Vte[t , «/], P^j (3.31) 



* P (t,^)-e ii [^(t / )] 



Pp 



> Vie [to, t/], pePy (3.32) 



which are recognized as conditions describing a DUE with elastic demand. 
3.3 The Fixed-Point Algorithm for Computing E-DUE 



□ 



Naturally Theorem 3.3 suggests the algorithm 



X k+i 



P Ao 



X k - aF{X h 



That is, at the k-th. iteration, we need to solve the following 



X 



k+l 



argrnin 
5 X 2 



X k - aF{X k 



X 



:IGA, 



(3.33) 
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which is recognized as a linear-quadratic optimal control problem 

^ +1 - (h k+ \ y k+1 (t f )) = a, 



X h 



dt 



s^ 11 f S \ E E [ h i^ - a *p^ hk ) - M*) 

+ J E (»6(*/) -«®5 [»*(*/)] -»«(*/)) a ( 3 - 34 ) 



yew 



subject to 



dyi 



dt 



E M*) V(i,j)GW 

y^(to) = V(i,j)eW 

h > 



(3.35) 

(3.36) 
(3.37) 



As we have already explained in the proof of Theorem 3.3 the optimality conditions for system 

(3.38) 



(3.33)-(3.38) are: 



A 



>.i 



y^(tf) + *%{y k (tf)] + yfr\t f ) v(», j) g w 



o < /^(i) 



[dH k 



arg 



= aig{-(^(t)-a* p (t, ^)-/ lp (t))-4(t / ) + aer.[/(t / )]+4+ 1 (t / ) = o} (3.39) 
In other words, given (i, j) G W and p G T^j, h k+1 {-) is determined as 



where 



K + \t) - ^ fc - «*,(*, h k ) + yj(t/) - «e-[/x*/)] - y-+\t f ) 

h k (t) - a%(t, h k ) + Q* + a9[Q fe ] - Q*+ 



Qij = !/§(*/) = E [ f h k p (t)dt V(i,j)€W 



(3.40) 



pePy 



Notice that for all (i, j) G W, Q 4 + must satisfy 



'7 



qT = E rV'w* 



pe7»« ^ 



£ / h k (t) - a^ p {t, h k ) + Q% + aQij [Q k ] - Q 



k+i 

ij 



dt 



(3.41) 



Equation (3.41) is viewed as a nonlinear equation with the unknown variable Q-J~ . The next 



lemma ensures that such equation always admits a unique solution. 
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Lemma 3.4. Given any ht{-) > 0, * p (-, h k ) > 0, Q k f > 0, &ij[Q k ] > and a > 0, there 



exists a unique Q- • G 3ft+ t/iat satisfies the following equation 






t f . 



fo 



/ij(*) - a¥ p (i, /i fc ) + gj + a0 o - [Q k ] - Qfr 1 ] dt - Q k + l = (3.42) 



Proof. We consider two cases. 

(i) h k p {t) - a^pitj^) + Q% + a@ i:i [Q k ] < V«(i) G [to, t/], Vp G ^. Then clearly g*+ J = 



is a solution to (3.42) 



(ii) h k (t) - a%(t, h k ) + Q k - + a0ij[Q fc ] > for some p G Vij and for t G B C [i , t/] where 



B is a set with positive measure. We call the left hand side of (3.42) f(Q;f ), which is a 



"'./ 



continuous function of Q- ■ . According to our hypothesis, the following hold 



/(0) > 0, 



/(Q 



Zc+In 



< if Q*? x » 1 



fc+i 



'./ 



where the notation 3> 1 means 'significantly large'. Therefore, by the Intermediate Value 
Theorem, there must exist at least one value of Q,- such that f{Q i t ) vanishes. 



'« 



The uniqueness of such a solution follows by observing that /(•), as a function of Q^ , is 



strictly decreasing. 



□ 



Lemma (3.4) suggests that one can conduct a simple root-search procedure to find the 
value of Q- ■ that satisfies (3.41 ). With preceding discussion, we can now state the algorithm 



for computing DUE with elastic demand. 
The Fixed-Point Algorithm for E-DUE 



Step 0. Initialization. Identify an initial feasible solution (h°, Q°) G Aq and set iteration 



counter k = 0. 
Step 1. Find the effective path delay. Solve the dynamic network loading problem with path 



flows given by h , and obtain the effective path delays \£p(-, h ), Vp G V . 



Step 2. Update the trip matrix. For each (i, j) G W, solve the following equation for Q^ , 



using root-search algorithms. 



Y, / [>»}(*) - «* P (t, &*) + Qfe + «e 4j [Q fc ] - Q*+ x j dt - Q^ +1 = 



pe7»« - 7 * 



Step 3. Update path flows. For each (i, j) G W and p G 7-V/, update the path flows as 
h k+1 (t) = \h k Jt) - a%(t, h k ) + Q% + a@[Q k ] - Q k + l ] Vt G [t , t f ] 



Step 4. Check for convergence. Terminate the algorithm with output X* ~ (h , Q ) if 



\(h k+l ,Q k+l )-(h k , Q k )\\ 
\\(h k , Q k )\\ E 

13 



< e 



where e G 3?++ is a prescribed constant, and 
\\(h,Q)\\ E = 



"£ I*' h 2 p (t) dt + Y, Ql) 
per Jto (i,j)eW J 



1/2 



is the induced norm on the product space E = (L 2 [to, t/]) x 3ft' ' . Otherwise, set k = k + 1 
and repeat Step 1 through Step 4. 



4 Bounded Rationality Dynamic User Equilibrium 

The notion of bounded rationality (BR) is a relaxation of the dynamic extension of Wardrop's 



first principle (Wardrop, 1952). The Wardropian principle enforces equal travel costs among all 



utilized route and departure time choices between an origin-destination pair. The bounded ra- 
tionality dynamic user equilibrium (BR-DUE), on the other hand, requires that the experience 
travel cost, including early and late arrival penalties, are within the interval [vij(h), Vij(h)+£ij]. 
Here Vij(h) is the essential infimum of the effective delays between origin-destination pair 
(i, j). Notice that such infimum depends on the complete path flow vector h. g^,- G K+ is a 
prescribed constant describing acceptable differences in travel costs experienced by travelers 
between origin-destination pair (i, j). 

In this paper, we present both variational inequality (VI) and differential variational in- 
equality (DVI) formulations of the continuous-time BR-DUE problem, when both route choice 
and departure time choice are employed. These formulations are original in the existing lit- 
erature. Moreover, the VI formulation provides a natural framework in which existence of 



BR-DUE can be analyzed by means of Brouwer's fixed point theorem; see Browder (1968) 
and Han et al. (2013c) for more detailed discussionF] The DVI formulation, on the other 
hand, allows the still emerging theory of differential variational inequalities to be employed 
for the analysis and computation of solutions of the BR-DUE problems. 

Throughout this section, the travel demand Qij between each origin-destination pair (i, j) 
is assumed to be a fixed constant. We begin with articulating the notion of bounded rationality 
in a measure-theoretic context. Recall from Section [2] that the effective delay operator \& maps 
each vector of path flows h G Aq to a vector of effective path delays, that is, 



*0, h) 



tf 



P\ ! 



h): pGV), %(;h): [t Q , t f ] -> K ++ Vp G V 



The set of feasible path flows, as we recall, is 



An 



h> 0: 



E 



f tf h p (t)dt = Q i3 V(i, j)ew\ c (L* ft,, t f ]) 



\v\ 



(4.43) 



The following definition of simultaneous route-and-departure choice BR-DUE is an extension 
of the static BRUE studied in the current literature: 

Definition 4.1. (BR-DUE) Given e = (e^ : (i, j) G W) G 3l + , a vector of departure 
rates h* £ Aq is a bounded rationality dynamic user equilibrium associated with the vector of 
tolerances e if for all (i, j) G W, 



h*Jt) > 0, peVi 



* p (t, h*) G [vijQf), Vij{h*) + Eij] V„(t) G [to, tf] (4.44) 



1 When we refer to the existence of BR-DUE, we are mainly concerned with solutions other than the dynamic 



user equilibrium in the sense of Friesz et al. ( 1993 1 



14 



where Vij(h*) is the essential infimum of the effective path delays between origin- destination 
pair (i, j). We denote this equilibrium by BR-DUE^, e, Aq, [to, tf]). 



4.1 The Variational Inequality Formulation of BR-DUE 

Essential to the variational inequality formulation of BR-DUE is the following operator: 

* £ : (L\[t ,t f ]f l -»• (L 2 + [t ,t f ]f l , h ^ (%(;h): p G V) (4.45) 

where 

$ 6 p (t, h) = max{* p (t, h), Vij(h) + eij} VpeTij (4.46) 

Given any vector of path flows h G Ao, the effective path delays Vl/p(£, h), p G V and the mini- 
mum effective delays Vij, (i, j) G W are uniquely determined by the dynamic network loading 
(DNL) procedure. Thus the operator & £ stated above is well-defined and strictly positive. We 
indicate the dependence of such an operator on the prescribed tolerance e by a superscript. 

The next theorem casts the BR-DUE^, e, Ao, [to, tf]) problem as an infinite-dimensional 
variational inequality problem. 

Theorem 4.2. (BR-DUE equivalent to a variational inequality) Given e G K^_ , let 
«K(-, h) : [to, tf] —7- K++ be measurable and strictly positive for allp G V and all h G Ao- Then 
a vector h* G Ao represents a BR-DUE^, e, Ao, to,tf) if and only if it solves the following 
variational inequality 



find h* G Ao such that 
Y I S %{t, h*){h p - h* p ) dt > 



P ev Jt ° 

V/iG A 



> VI($ £ , A , [to, t f ]) (4.47) 



Proof, (i) [Necessity] Let h* G A correspond to a BR-DUE solution, then for any h G A and 
any (i, j) G W, 



E E / f %& h*)h p dt > y E / ' 4^*) ^ dt 



where 
and 



= E 4-( /i *) E / h p (t)dt = y v-u h *)Qij 

(4.48) 

/4j(fr*) = min {fj, E p (h*) : p G 7^} 
^(/i*) = essinf{$£(t, /i*) : i G [t , tf]} 
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It is easy to deduce that /if(/i*) = Vij(h*) + S{j. On the other hand, with (4.44) and (4.46) 
we perform the following calculation: 



E E f f %(t,h*)h* p dt = Yl E /^ {vii(h*) + eij)Kdt 

= E (Mh*)+eij)Q*j = E tfj( h *)Qii ( 4 - 49 ) 



(i,j)ew P ev 1J Jt ° 



(i,i)evv 



(M)eW 



In view of (4.48) and (4.49), we have 
ft 



E E [ * $p(t,h*)h* p dt < E E [ f %(t,h*)h p dt V/iEA 



(i,j)&wpeVi; Jto 



which is recognized as the variational inequality (14.47) 



(ii) [Sufficiency] Let /i* E A be a solution of U/(<I> e , A , [to, tf]). Given that $ e : (L+[£ , t/]) |P| 

(L^_[to, £/]) is measurable and strictly positive, we employ the same proof of Theorem 2 in 
|Friesz et al.| ([1993) to demonstrate that h* satisfies 



h* p (t) > 0, peVij => $|(t, /»*) = ^-(/i*) 
That is, 



V„(t) E [to, t f ], V(i, j) 6 W (4.50) 



/i!(t) > 0, p E 7\ 



%(/i*) < *£(*, h*) < v l3 {h*)+e V„(t) E [t , t f ), V(i, j) E W 



This concludes that /i* solves the BR-DUE problem. 



D 



Let us illustrate the meanings of VI (\P, A , [to, tf]) (DUE) and Vl(® £ , A , [to, tf]) (BR- 
DUE) graphically. The following discussion is mathematically less rigorous but should provide 
an effective visualization of the problems of interest. Let us consider a network with just one 
link and one origin-destination pair. Figure [T] depicts a possible scenario of a DUE path flow 
h*(-) and the associated effective path delay ^ p (-, h*). It is intuitively clear that in order to 
minimize the following quantity over all h p 



* p (t, h*)h p dt, 



(4.51) 



to 



the non-zero portion of h*(t) should be located at the "flat bottom" of the effective delay 
curve. Now let us turn to BR-DUE, which requires that whenever h*{t) is non-zero, 9 p (t, h*) 
must be within an indifference band with width e. Clearly the upper boundary of such band 
is Vij{h*) +£. Thus being a BR-DUE path flow requires that the non-zero portion must reside 
within the time interval [a, b] (see Figure^); i.e. the BR-DUE path flow maximizes 



to 



$p(t, h*)h p dt 



(4.52) 



Theorem 4.2 presents a generic variational inequality that is equivalent to the BR-DUE 



problem. In comparison with the VI formulation VIC&, Aq, [to, tf]) of SRDC DUE (Friesz et 
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%(tX) 



\V \ 




Figure 1: A possible scenario of a dynamic user equilibrium path flow h*(-) and the associated 
effective path delay ^ p (-, h*). In order to minimize the quantity (4.51) over all h p , the equi- 



librium path flow h* must be located at the "flat bottom" of the effective delay curve. Such 
observation also coincides with the definition of DUE. 



h;© 



%(t,h*) 




Figure 2: A possible scenario of a BR dynamic user equilibrium path flow h*(-) and the 
associated effective path delay ^ p (-, h*). To be a BR-DUE, h* must vanish outside the time 
interval [a, b], which implies that h* is located at the "flat bottom" of the curve formed by 
the graph of ^fp(-, h*) and the horizontal line Vij + e. Such a curve is precisely the graph of 
3> £ (-, h*). Therefore, the BR-DUE is an analogy of the DUE when the operator \& is replaced 
with <£ £ . 



al. 



1993), the proposed variational inequality Vl(& 6 , Aq, [to, if]) has a different principal op- 



erator <i> e , which encapsulates the dynamic network loading procedure as well as the tolerance 
vector e. However, these two types of Vis do share a similar mathematical structure which 



allows known methodologies regarding SRDC DUE, conveyed by a series of papers (Friesz et 



al., 1993, 2001, 2013 2011; Friesz and Mookherjee 2006 Han et al., 2013c), to be transferred 



to BR-DUE problems. In fact, we are able to establish the differential variational inequal- 
ity (DVI) formulation and the fixed-point algorithm for BR-DUE, as is demonstrated in the 
remainder of this section. 
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4.2 The Differential Variational Inequality Formulation of BR-DUE 



Let us introduce the function yij(-) : [to, tf] — > [0, Qij] for each origin-destination pair (i, j); 
Dij{t) represents the total traffic volume that have departed from origin i with the intent of 



reaching destination j. Note that the set of feasible path flows Ao defined in (4.43) can be 
re-stated as 



A, 



h> 



dy 



ij 



dt 



E M*)> unto) = 0, yij(t f ) = Qij V(»,j)eW 



(4.53) 



As a result, we have the following useful theorem. 



Theorem 4.3. (BR-DUE equivalent to a differential variational inequality) Given 
e G W + , assume $ e (-, h) is measurable and strictly positive for all p G V and all h G A2. A 
vector of departure rates (path flows) h* G A2 is a bounded rationality dynamic user equilibrium 
if and only if h* solves the following differential variational inequality 



find h* G A2 such that 

E [ f %(t,h*)(h p -h* p )dt > 



P&V 



V/lGAa 



► DVI{$ £ , A 2 , [t , */]) 



(4.54) 



Proof, (i) [Sufficiency] Given that h* G A2 is a solution of the differential variational inequality 
DVl(& £ , A2, [to, t/]), we derive the following obvious relation: 

E I ' *k*, >»>;(*) * < E / ' $ p(*> ^)M*) df 

for all /i G A2. In other words, /i* is the solution of the optimal control problem: 



[tf 
min J(h) = E E / %(t, h*) h p (t) dt + ^ My[Qii- Jfcj(*/)] 



(4.55) 



(iJ)eWpeVij to 



(i,i)ew 



subject to 



eft 



yy(t) = ^ ^(t) V(», j) G W 

Wi(to) = V(t, j) G W 
/i > 



where each //jj is the dual variable for the terminal conditions on yij( 
Hamiltonian for problem (4.55|)-p.58[) is 



(4.56) 

(4.57) 
(4.58) 

(i, j) G W. The 



H = E E w»*)^+ E ^ E ^ 

(j,i)eWpeP<> (i,i)ew pe7\,- 



where the adjoint dynamics are: 



d_ 
dt 



hjit) 



dH 

dyij 



V(», j) G W, t€ [to, tf] 
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(4.59) 



the transversality conditions read: 

. f s ®Yl(i,j)eW Vij [Qij -yij(tf)\ 

Taken together, (4.59) and ( 4.60| ) imply that 

Xij(t) = - nij V(t, j) g W, Vt € [to, */] 

According to the minimum principle, we have 



V(i, j) EW, t G [to, t/] (4.60) 



(4.61) 



/l* 



arg min H such that — /i < 



for which the Kuhn- Tucker conditions are 

%{t, h*) + \j = Pp > Vp G 7^-, V(t, j) G W, V„(t) G [t , t/] 

Pp (t) h* p {t) = Vp G P 4j , V(», j) G W, V„(t) G [t , t f ] 

P P (t) > VpeVij, V(», j)GW, V„(t) g [t , t/] 



(4.62) 
(4.63) 
(4.64) 



where p p are dual variables for the non-negativity constraints. From ( |4.62[ ) and (4.63) we have 
that for all (i, j) G >V, 



/i;(t) > 0, p G Py => %{t, h*) = - Xij = pij V„(t) G [t , t f ] 
$p(t, h*) = Pp-Xij = pp + pij > mj M v {t) G [t , tf] 

from which we conclude that 

IHj = essinf{^(t, h*) : t£ [t , t/]} V(i, j) G W 



(4.65) 
(4.66) 



In view of the definition of <K (4.46), we re-state ( 4.65[ ) as 

h* p (t) > 0, p G 7> 4j => *(*, /i*) < %(&*) + £« V„(t) G [t , t 7 ], V(i, j) G W 

which is recognized as conditions describing a dynamic user equilibrium with bounded ratio- 
nality, where Vij(h*) is the essential infhnum of effective path delays between origin-destination 
(hO)- 

(ii) [Necessity] For each origin-destination pair (i, j) G W, define pfj(h) to be the essential 
infimum of $£(•, h) between (i, j). It is then easy to verify that nh(h) = Vij(h) + £y for all 
(i, j) G W and for all /i G A 2 . 

For a BR-DUE solution /i* we have for all (i, j) G W that 



which translates to 

/»;(<) > o, P g Vij 



Vij(h*) < ¥ p (t,h*) < Vij(h*) + Si 



%(t, h*) = pfjih* 



V(i, j) G W 
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dt 



(i,j)ew P £P tj 



Therefore, given any h £ A2, we have 

E E f f %(t,h*)h* p {t)dt = y, E f^-wv*) 

E /W E /*'*;(*)* 

E ^W)-Qv 

(i,i)evv 

E ^#*) E f f h P {t)dt 



(i,j)6VV 






^ E E [ f %(t,h*)h p (t)dt 



Therefore, /i* G A 2 is a solution of the DVI (4.54). 



□ 



4.3 The Fixed-point reformulation and algorithm 

Not surprisingly, the BR-DUE also admits a fixed-point-problem reformulation provided the 
differential variational inequality established in Theorem |4.3| 



Theorem 4.4. (Fixed point problem equivalent to BR-DUE) Assume that ^ p (-, h) is 
measurable and strictly positive for all p £ V and h £ A2. Then the fixed point problem 



h* = P A2 [h* -a<S> £ (t, h*) 



(4.67) 



is equivalent to DVl(& e , A2, [to, £/]) and hence to BR-DUE^, e, Ao, [to, £/]), where P\ 2 
is the minimum norm projection onto A2 and a > is a fixed constant. 



The proof of Theorem 4.4 follows closely that of the fixed-point reformulation of SRDC DUE 



provided by Friesz et al. (2011). The only difference is that the effective delay operator ^ 



therein is now replaced with the newly defined <£ £ . Thus we omit the proof from this paper. 
The purpose of invoking the DVI formulation and the fixed point problem formulation is to 



devise an iterative computational scheme. Namely, Theorem 4.4 suggests the algorithm 



h 



fc+i 



Pa 2 



h k - a$ £ (t, h k 



That is, at the k— th iteration, one needs to solve the following linear-quadratic optimal control 
problem 

1 



h 



k+l 



argrnin 
h I 2 



h k - a$ e (t, h k ) - h 



h £ A- 



That is, we seek the solution of the optimal control problem 

2 L 



mm J k {h) = Yl IH j [Qij-Vij(tf)]+ E E/ ^[h k p -a%(t,h)-h p \ (4.68) 



(i,j)ew 



{i,j)&WpeV i: 
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subject to: 



dyij 
dt 



Vij(t 



h > 



£M*) V(i,j)eW 

V(i,j)eW 



(4.69) 

(4.70) 
(4.71) 



Finding dual variables associated with terminal time demand constraints turns out to be 
relatively easy. Note that the relevant Hamiltonian for (4.68), (|4.69[), (|4.70[) and (4.71) is: 



» k = \ E T,[K 



(i,j)eWp£Vij 



' a%(t,h k )-h p 



+ 2^ ^J 2^ h p 



(4.72) 



where each A,,- is an adjoint variable obeying 



dXij 
~dT 



, ,dH k 
(-1 ) lr - = 

d iMj [Qij ~ Vij (t)] 



V(t, j) G W 

= -^j V(i,j)eW 



9yij(tf) 

From the above we determine that 

Xij(t) = - fiij Vt G [to, tf], (i, j) G W 
The minimum principle implies for any p G V, 



h k+1 (t) 


= arg 


I un 

I dh p 


'hus 


we 


obtain 








h k+1 


(t) = 



z\j§~ = 0J = arg{[(^(t)-a^^,/ l ^_/ lp (t))](_l)_ /Xu . = } 



h;{t)-a<S> p [t,h k \ +fiij 



V(i,j)eW, PGP, 



(4.73) 



Notice that the following flow conservation constraint applies here. 



£ f 1 h k+1 (t)dt = Qi 



V(i,j)GW 



pern ° 
Consequently the dual variable v^ must satisfy 

nt f 



E 



peVij Jt ° 



h*(t)-a$Ut,h 



Vij 



dt 



Qi 



V(t, j) G W 



(4.74) 



Recall that each fiij is time invariant and notice that the equations of (4.74) are uncoupled; 
furthermore, the left hand side of (4.74) is strictly increasing in mj, we see that simple root 
search will find the value of \iij satisfying (4.74). Once the fJ,ij' s are determined, the new 
vector of path flows h +1 for the next iteration may be computed from (4.73). As a result, we 
have the following algorithm. 

Fixed-point Algorithm for BR-DUE 
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Step 0. Initialization. Identify an initial feasible solution h° G A2 and set iteration counter 



k = 



Step 1. Find <3? £ . Solve the dynamic network loading subproblem with path flows given by 



h k , and obtain the effective path delays ^>(-, h k ). Then define 

$ £ p (t, h k ) = max{* p (i, h k ), %(/i fc ) + 6^} 
where v- L j{h ) is the minimum effective delay within origin-destination pair (i, j) 



using root-search algorithms. 



Step 2. Find the dual variable. For each (i, j) G W, solve the following equation for fiij, 



E 



to 



h k (t)-a$ e p (t,h k ) +Hij 



dt 



Step 3. Update path flows. For each (i, j) G W and p G Vij, update the path flow as 



,fc+i 



K +l (t) 



h;(t)-a%[t,h K ) +ntj 



Vt G [to, t f ] 



Step 4. Check for convergence. Terminate algorithm with output h* ~ h if 

\\ h k+i_ h k\\ 



\l 2 



\\h k \\ L2 
where e G K ++ is a prescribed constant, and 



< e 



h\\ L 2 = Yl [ f h 2 p {t)dt 
P ev •' to 



is the standard L 2 -norm on (L 2 [£o> tf]) . Otherwise, set k = k + 1 and repeat Step 1 through 
Step 4. 

5 Existence of solutions 

The most obvious approach to establishing existence for any of the user equilibrium models 
presented above and in the existing literature is to convert the problem to an equivalent fixed 
point problem and then apply Brouwer's fixed-point existence theorem. One statement of 



Brouwer's theorem appears as Theorem 2 of Browder (1968): 



Theorem 5.1. (Browder, 1968) Let K be a compact convex subset of the locally convex 



topological vector space E, T a continuous (single-valued) mapping of K into E* . Then there 
exits uq in K such that 

(T(u ), uq-u) > 

VuG K. 



Proof. See Browder (1968). 



□ 



Approaches based on Brouwer's theorem require the set of feasible path flows (departure 
rates) under consideration to be compact and convex in a Hilbert space, and typically involve 



an a priori bound on all the path flows (Zhu and Marcotte, 2000). 
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5.1 Existence of E-DUE 



Analysis regarding the existence of E-DUE is most facilitated by expressing it as a variational 
inequality: 



find X* £ A such that 



(F(X*),X-X*) E > 



Vie A 



>vi(F,e,A,to,t f 



(5.75) 



where X = (h, Q), A is defined in (3.10), and J-(-) is defined in (3.19). The inner product 



(•, -) E on the extended s pace E = (L 2 [t , t f ]y ' x Kl w l is defined in ( |3.17| ). Unlike the 



3.16 



5.75 



SRDC DUE introduced by Friesz et al.| ( 1993"| ), the DVI representation of E-DUE expressed by 

has both infinite-dimensional and finite-dimensional terms. One advantage of invoking 

is a more compact mathematical expression which absorbed these inhomogeneous terms. 

herefore, in order to apply Theorem 5.1, one only has to investigate the operator J-(-) and 

the subset A. 

Nevertheless, the aforementioned approach faces several theoretical difficulties: (1) the set 
A is not compact in E; and (2) a direct topological argument requires a priori bounds for the 
path flows, where those bounds do not arise from any behavioral argument or theory. These 



two difficulties are overcome by Han et al. (2013c) in the derivation of existence result for SRDC 
DUE with fixed demand. In that paper, the first difficulty is handled by successive finite- 
dimensional approximations of the set Ao; while the second one is overcome by introducing 
minor assumptions on the travlers' disutility function. 



Using techniques introduced by |Han et al. (2013c), we may show that the variational 
inequality (5.75) has a solution, a detailed proof is provided by Han et al. (2013d). 



5.2 Existence of BR-DUE 

We are mainly concerned with the existence of BR-DUE solutions that are not DUE solutions, 
since the latter is trivially a BR-DUE solution. Before the variational inequality formulation 
is introduced in this paper, very few analytical tools are available for studying existence of 



BR-DUE. Szeto and Lo (2006) propose to analyze the existence of route-choice BR-DUE 



by introducing a gap function, which measures the differences in travel costs among various 
route choices. However, that approach does not provide conditions under which a solution is 
guaranteed to exist. 



Let us recall the VI formulation of SRDC BR-DUE established in Theorem 4.2 According 



to Theorem 5.6 of Han et al.l (|2013c), such a VI has a solution provided that the operator 

\T>\ VU\ — 

<3? e : (L 2 [£o, £/]) — > (L 2 [tQ, tjfj is continuous, together with some other minor regularity 
assumptions. In general, <& £ can be continuous even if the effective delay operator \I> is not. 
Therefore, conditions that guarantee the existence of BR-DUE should be weaker than those 
for DUE. 

In order to say something about the continuity of <3? E , one needs to investigate the under- 
lying dynamic network loading (DNL) model. In particular, we have the following existence 



theorem when the DNL procedure employs the Vickrey model (Vickrey 
2013a|b| ) 



1969, Han et al. 
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Theorem 5.2. (Existence of BR-DUE) Let the effective path delay be defined as 

%(t, h) = D p (t, h) + f(t + D p (t, h) - T A ) Vp£V 

where h € Aq, T^ is a fixed constant indicating the target arrival time. Assume that /(•) is 
continuous on [to, tf] and satisfies 

/(fe)-/(ti) > A(t2-ti) Vt < h < t 2 < t f (5.76) 

for some A > — 1. Furthermore, assume that the dynamic network loading procedure em- 
ploys the Vickrey model to describe flow dynamics. Then, For any fixed e £ ffl ', BR- 
DUE^, e, Aq, to, tA as in Definition 



has a solution. 



3> £ 



Av\ 



Proof By Theorem 5.6 of Han et al. (2013c), it suffices to show that the operator 

(L 2 [to,t f ]) m 

VpeVij, V(i,j)eW 



{L 2 [t , t f ] 
$ £ p (t, h) = max{$ p (t, h), %(/i) +%} 



(5.77) 



is continuous. Let {h n } C Aq be an arbitrary strongly converging sequence with limit h. By 



Theorem 5.3 of Han et al. (2013c), the corresponding effective path delays ^(t, h n ) converges to 



ty(t, h) uniformly for all t 6 [to, tf]. Thus the function $ e (-, h n ) defined in ( |5.77 ) converges to 
$ £ (-, h) uniformly and also in the L 2 -norm. This shows that the operator <J> £ is continuous. □ 



Szeto (2003) identifies certain dynamic network loading procedures that do not admit a 



continuous effective delay operator. With these DNL models, a DUE in the sense of Friesz 
et al. (1993) may not exist. However, if one can show that the corresponding operator <l> e 



is continuous, then the existence of a BR-DUE can be guaranteed by invoking the analytical 
framework proposed in this section. Showing the continuity of <J> e when the corresponding 
effective delay operator is not continuous is beyond the scope of this paper but will remain an 
interesting research direction for traffic modelers. 

6 Convergence of the fixed-point algorithms 



Both fixed-point algorithms presented in Section [3. 3| and Section [4. 3| are particular instances 
of the abstract algorithm 



v fc+i 



G(z k 



(6.78) 



for solving the fixed point problem 



G{z*) 



z* G V 



where V is a proper Banach space. If the abstract operator G(-) satisfies the contracting 
property, i.e. 

\\G(x) - G{y)\\ < 6\\x-y\\ Vx, y G V (6.79) 



for some < 5 < 1, then the sequence generated by (6.78) converges to the fixed point z* 



according to the Contracting Mapping Theorem (Rudin, 2006). However, the contracting 



property (|6.79|) is too restrictive in application, thus we need a relaxation known as non- 

G(y)\\ < \\x — y\\ \/x,y£V 



expansiveness, that is, 



G(x) 



(6.80) 



The following theorem, taken from Friesz et al. (2011 ), is particularly useful to establish strong 



convergence of algorithms in the presence of non-expansive G(- 



21 



Theorem 6.1. Consider the abstract mapping G : V — > V, where V is a Banach space. The 
sequence generated by 

z k+1 = (3 k z° + (l-(3 k )G(z k ) (6.81) 

converges strongly to a fixed point z* of G(-), that is, 



lim z k 

k— »oo 



G(z*) 



when G(-) is non- expansive, z° £ V is an arbitrary point and 
(i) fa 6 [0, 1] 
(ii) limfc^oo f3 k = 

(Hi) Efclo h = oo 

(iv) lim fc ^ 00 (/3 fc - P k _i)(f3 k y l = 0. 



Proof. The reader is referred to Friesz et al. (2011) for a proof of strong convergence. 



□ 



We call (6.81) the modified fixed point algorithm. By virtue of Theorem 6.1, in order to 



analyze convergence of the proposed fixed point algorithms, it suffices for us to check that the 
operator G(-) is non-expansive. 

6.1 The E-DUE case 



Let us recall the fixed-point iteration scheme from Section 3.3 



X K 



P ~Ao 



X K - aT{X h 



k = 0, 1, ... 



where P-r [■} is the minimum-norm projection onto the subset Aq, X = (h, Q) E Aq, and 



F{X) = T{(h,Q)) = (*(-, h), e~[Q]) 



where G [Q] = —Q[Q]. It is well-known that projection operators of the type P-r [•] are always 
non-expansive. We have the following convergence result. 

Theorem 6.2. (Strong convergence of the modified fixed point algorithm for E- 

DUE) Assume that the effective path delay operator ^ and the inverse demand function Q 
satisfy the Lipschitz condition, that is, 



|tf(-, h 1 )-^, h 2 )\\ L2 < L^h 1 - h 2 \\ L2 \fh\h 2 e(L 2 + [t ,t f ]) lvl 
||e[Qi]-e[Q 2 ]|| < L 2 \\Qi-Q 2 \\ \/Q 1 ,Q 2 eR 1 ™ 1 



(6.82) 
(6.83) 



, 1^1 



for some L\ > 0, L 2 > 0; and ^ is strongly monotone on (L^fto, tf]) with constant K\ > 0, 
i.e., 



<*(•, h l )-V(-, h 2 ), h 1 -^) > IdW^-h 2 " 2 



lL2 Vh 1 , h 2 G (L 2 + [t , tf}) 

In addition, assume that the inverse demand function satisfies 



\v\ 



(@[Qi] - 6[Q 2 ], Qi - Q 2 > < - K 2 \\Q Y - Q 2 \\ 2 VQi,Q 2 £^ 1 
for some K 2 > 0. Then the modified fixed point algorithm converges strongly. 



(6.84) 



25 



Proof. According to the preceding discussion, we only need to show that the map defined as 



G(X h 



P A 



X k - aF(X l 



is non expansive. Since the minimum norm projection is non-expansive, we have 
rk+1 )-G(X k ) ' 

E 

-k i „ T~rvk\ 



< 



G(X k+1 ) - G(X' 
fc+l „,Tfvk+l' 



X K+l - aT(X K+i ) - X* + aT(X 

2 



X 



fc+l 



X 



fc+l 



X k 



x k 



2a(F{X k+l )- F{X k ), X k+l -X k ) +a 2 F{X k+l ) - F{X k ) 

E \ IE 

2 



2al (*(-, h k+l ) - *(-, h k ), h k+l - h k ) + (®-[Q k+l ] - Q-[Q k }, Q k+1 - Q k ) J 



+ a 



*(■, /i fc+i )-^(-, /i fe ) 



L 2 



+ 



< 



X 



fc+l 



X A 



2a ifi 



,fc+i 



fc A 



-[Q^+i]_e-[Q fc ] 

Q k+1 - ( 



+ X 2 



+ a 2 \L\ 



,fc+i 



+ L| 







fc+l 



< 



x 



fc+i 



x A 



2amin \Ki, K2] 



I? 



Q l 



vk+1 -y-k 



+ a max {L 1 , L 2 } 



V'fe+l -y-k 



X 



fc+l 



= (l-2amin{ETi, ET 2 } +a 2 max{Li, L2}) 
Thus, to ensure non-expansiveness of G(-), it suffices to let 

(l -2amin{Ki, K 2 } + a 2 max{Li, L2}) < 1 = 



X" 



a < 



2mm{K 1 ,K 2 } 
max{L 2 , L|} 



D 



Remark 6.3. The convergence issue for computing dynamic user equilibrium has been a chal- 
lenge for decades. Difficulties arise mainly from the non-monotone nature of the path delay 
operator or the effective path delay operator. Very few dynamic network loading models are 
reported in the literature that yield a strongly monotone delay operator, except the ones pro- 
posed by Mounce \200(fy and Perakis and Roels (200(fy , which rely on somewhat restrictive 
assumptions. The contribution claimed by Theorem \6.2\ should be interpreted as a generaliza- 
tion of the convergence proof for the fixed demand case to an elastic demand case, rather than 
solving the problem of monotonicity of the delay operator, which of course, is quite desirable 
and remains a focus of future research. 



The assumption (6.84) made on can be illustrated as follows. At first glance it suggests 



that —0 should be strongly monotone with constant Ki > 0. If we can further assume that 
for each (i, j) G W, Vij = ®ij[Qij], i-e. the inverse demands of different origin-destination 
pairs are decoupled, which is reasonably true, then (6.84) reduces to the following: 






< 



K 



1:1 



VQij > Qn, V(i, j) e w 
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or simply 



G'iiiQ 



ijr'LVJJJ 



< 



Ki. 



V(i, j) G W 



if ®ij is differentiable, where Kij G 3?++ is a constant depending on origin-destination pair 
(i, j) G W. Clearly, such assumption holds for any afnne or exponential functions. 



6.2 The BR-DUE case 

Regarding simultaneous route-and-departure choice dynamic user equilibrium with bounded 
rationality, we have the following convergence theorem for the fixed-point algorithm. 

Theorem 6.4. (Strong convergence of the modified fixed point algorithm for SRDC 

BR-DUE) Fix e G 3^' , assume that the operator $ E : (L 2 + [t , t/]) |P| :-)■ {L 2 + [t , t/]) |P| 
defined as 

$ £ p (t, h) = max {*(*, fc), Wii (/i) + Sij) Vt G [to, t f ], p£ P ijt (i, j) G W 

satisfies the Lipschitz condition with constant L > 0; and <J> e is strongly monotone on Ao mi/i 
constant K. Then the modified fixed point algorithm converges strongly. 



The proof of Theorem 6.4 follows closely that of Theorem 6.2 , and is omitted in this paper. 



7 Numerical studies 

In this section, we will test the numerical performances of the proposed fixed-point algorithms 
for computing E-DUE and BR-DUE. Regarding the dynamic network loading (DNL) model, 



we employ the LWR-Lax model proposed by Friesz et al. (2013). Nonetheless, the fixed 



point algorithms do not rely on any specific forms or assumptions of the link dynamics, flow 
propagation constraints and delay models; in other words, they can be equally applied to all 
DNL models. 

All of the reported computations of the fixed-point algorithms are performed on a 2.7 GHz 
4 GB of RAM computer. Every attempt was made to avoid the use of numerical tricks that 
could not be equally applied to all models. 

7.1 DUE with elastic demand 

We conduct numerical study of the fixed-point algorithm for E-DUE on two networks: a simple 
two-node, two-link network, and the 24-node, 76-arc Sioux Falls network. 

7.1.1 Simple network 

We begin with a simple network with just one origin-destination pair (1, 2) and two paths p\ 
and p2 , depicted in Figure [3j The link parameters are given in the table next to the network 
figure. 




Arc 


Jam density 
(vehicles / km) 


Free flow speed 
(km / 5min) 


Length 
(km) 


1 


800 


6 


4 


2 


600 


3 


4 



Figure 3: Topology and arc-specific parameters of the simple network. 
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For simplicity, we choose the inverse demand function as an affine function: 



V12 



6i 2 [Q 



12 



0.025Q12 + 100 



(7.85) 



The fixed-point algorithm took 33 iterations to declare convergence. The elastic demand, Q k 2 , 
at each iteration is shown in Figure |4j The average effective delay between origin-destination 
pair (1, 2) at each iteration is shown in Figure [51 We note that in each intermediate iteration, 
the effective delays experience by drivers are not the same over all departure time and routes. 
For this reason we can only show in Figure [5] the average effective delay computed as 



"12 



/// [*!(t, h k ) h k (t) + tt 2 (t, h k ) h k 2 (t)] dt 



The effective delays are gradually equalized when the algorithm converges. We can observe 
clearly from Figure UJ and Figure pi the converging patterns of Q k 2 an d v k 2 ■ Figure \6\ shows the 
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^000000000000000000000 
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Iteration Number 



Figure 4: Simple network: Elastic demand at Figure 5: Simple network: Average travel cost 
every iteration at every iteration 

final path flows and effective delays when the algorithm terminates. We are assured that not 
only do Q k 2 and v k 2 converge, but also the effective path delays become equilibrated. When 
this happens, the average effective delay also becomes the minimum travel delay. 





Figure 6: Simple network: Equilibrium path flows and effective path delays. 
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The relative distance between two consecutive iterates is denned as 

\\(h k +\Q k 



D 



rel 



,/•-!-: ni>-L\-(h k , Q k )\\ 2 E 



\\(h k , Q k )\\l 



(7.86) 



The fixed-point algorithms are terminated once D k el < e, where e is a fixed threshold. Figure 
7\ shows D k „, for each k, from which we confirm the convergence of the fixed-point algorithm. 
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Figure 7: Simple network: Relative distance D k el 



in percent) at each iteration. 



Finally, we need to check if v* 2 and Q\ 2 claimed by the fixed-point algorithm satisfy the 
prescribed inverse demand relation (7.85). Figure^shows the pair (Q k 2 , ^12) a ^ each iteration 
k. We observe that the generated sequence (Q k 2 , v k 2 ), k > 1 converges with the limit point 
on the straight line defined by the inverse demand function. 
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Figure 8: Simple network: Plot of the pair (Q k 2 , v k 2 ) at the /^-iteration. Notice that each v\ 2 
is the average travel cost of drivers. The solid line represents the prescribed inverse demand 
function. 
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7.1.2 The Sioux Falls network 

We test our proposed fixed-point algorithm on a larger network (Sioux Falls) depicted in 
Figure |9j This network consists of 24 nodes and 76 arcs. We consider 10 origin-destination 
(OD) pairs and 209 paths for our numerical study. The set of OD pairs reads 

W = {(1, 20), (3, 20), (4, 20), (5, 20), (9, 20), (12, 20), (11, 20), (10, 20), (2, 20), (6, 20)} 

Due to space limitation, we do not elaborate on our path set or the arc-specific parameters 
employed by the dynamic network loading procedure. 




Figure 9: The Sioux Falls network, with 10 origin-destinations pairs and 209 paths selected. 



We employ the following affine inverse demand function for all OD pairs: 



&ij[Q 



iji 



0.025Qij + 150 V(i, j) <= W 



(7.87) 



The fixed-point algorithm takes 36 iterations to claim convergence. The elastic demands 



between the 10 OD pairs at each iteration are shown in Figure 10, from which we observe 
converging patterns of these elastic demands. 

For each origin-destination pair (i, j) G W, we plot the pairs (Qh, vfA for k = 1, . . . ,36, 
where Q\a and v\ A are the travel demand and the average travel cost at the k -iteration, 



respectively. The results are summarized in Figure 11. It is clearly shown that each sequence 



(Qiji v ij)> (*> J) e W converges to some point on the line corresponding to the inverse demand 



function ( |7.87[ ). 

In order to demonstrate that when the algorithm terminates, the effective path delays 



within each OD pair are equalized, we show in Figure 12 the departure rates and the corre- 



sponding effective delays of two paths pso and piso- From Figure [12] we see that not only do 
Qf A and v\ A converge, the travel costs are also equilibrated and independent of departure time 



'*j 



tj 



selected by travelers. Finally, the relative distance D^ el defined in (7.86) is plotted against the 



iteration number in Figure 13 The relative distance is monotonically decreasing with each 



iteration, and satisfies the termination criterion (< 0.01) after 36 iterations. 
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Figure 10: Sioux Falls network: Elastic OD demands at each iteration. 
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Sioux Falls network: Convergence of the sequence {Q%, «£•) to the inverse demand 
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Figure 12: Sioux Falls network: Path departure rates and effective path delays. Left: Path 
50; right: Path 150. 
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Figure 13: Sioux Falls network: Relative distance D^ , (in percent) at each iteration. 

7.2 DUE with bounded rationality 

In the first half of this section, we will illustrate the solutions of the BR-DUE problem com- 
puted by the fixed-point algorithm on two networks: the simple network and the Sioux Falls 
network. In the second half, we will test the convergence of the algorithm on four different 
networks. 



7.2.1 Simple network 



Let us consider again the simple network illustrated in Figure [3| Assume there is a fixed 
demand Q12 = 3300 between origin-destination pair (1, 2). As termination criterion for the 
fixed-point algorithm, we stipulate that -D^f < 0.001 
between two consecutive iterates: 



where -D^f is the relative distance 



D 



ref 



\\h k+1 -h k \\ 



L^ 



\h k 



2 
\l 2 



We consider the following value of tolerance: e = 15. The fixed-point algorithm takes 11 
iterations to declare convergence. The resulting BR-DUE path flows and associated travel 



delays are shown in Figure 14. We can see that when the algorithm terminates, the path flows 



indeed satisfy the variational inequality (4.47) as they are located at the "flat bottom" of the 



curves of <I> £ . The second row of Figure [14| shows the actual effective delay. One can see that 
the travel costs experienced by all travelers are within the indifference band with width 15. 



7.2.2 The Sioux Falls network 

We consider the Sioux Falls network with the same set of origin-destination (OD) pairs and 
the same set of paths as Section 7.1.2 For simplicity, we assume the same tolerance e = 10 



for all OD pairs. The termination criterion is again D k e r < 0.001. 

The fixed-point algorithm takes 15 iterations. Figure 15 shows two path flows with their 
corresponding <I> e and \P. The first row shows flat bottom of the functions $25 an d ^200 1 
indicating that the path flows generated by the algorithm solves the variational inequality 



(4.47). The second row shows the actual effective path delays experienced by travelers, which 
are within the range of e = 10. 
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Figure 14: Simple network: BR-DUE when the tolerance e = 15. First row: Path flows and 
the corresponding 3> e (-). Second row: Path flows and the corresponding effective delays \l/(-) 

7.2.3 Convergence of the fixed-point algorithm for BR-DUE 

In this section, we test the convergence of the fixed-point algorithm with different values of e. 
For the comprehensiveness of our experiment, we consider two additional networks depicted 
in Figure [T6] , so that there are totally four networks: the two-arc, two-node simple network; 
the six-arc, five-node small network; the 19-arc, 13-node medium network; and the 76-arc, 
24-node Sioux Falls network. 
The termination criterion is, 



U ref 



\h 



fc+i 



h 



fell 2 

L 2 



\h k 



i2 



< e 



where we choose e = 0.001 for the first three networks and choose e = 0.01 for the Sioux Falls 
network. 

We test the proposed fixed-point algorithm with different values of the tolerance e. The 
results are summarized in Table [T] through Table [4j From these tables, we observe a dis- 
cernible improvement of convergence when the value of e becomes larger. This coincides 
with the intuition that when the indifference band is wider, a BR-DUE is more likely to be 
reached. However, we note that in practice, the choice of e should be based on a behavioral 
study of travelers and/or empirical calibration, rather than a numerical trick used to facilitate 
convergence. 



33 



•^ 10 



- 


h / : 


1 

Path 25 Departure Rate 



3 

c 12 
| 




n / 


J 
















■* 8 






of 








^ 


Path 200 Departure Rate 




- 


Q 


200 


| 



Time (minute) 



Time (minute) 



■■— 10- 

3. 




a 

"SOS 




Path 25 Departure Rale 
Path 25 Effective Delay 



I 4 " 

Q 



- Palh 200 Departure Rate 

- Palh 200 Effective Delay 



Time (minute) 



Time (minute) 



Figure 15: The Sioux Falls network: BR-DUE when the tolerance e = 10. First row: Path 
flows and corresponding <J> £ . Second row: Path flows and the actual travel costs. 





Figure 16: Left: The six-arc, five-node small network. Right: The 19-arc, 13-node medium 
network. 





e = 


e = 5 


£ = 10 


£ = 15 


e = 20 


e = 25 


Iteration 


19 


16 


15 


11 


9 


9 


Time (s) 


3.09 


2.72 


2.49 


1.91 


1.58 


1.61 



Table 1: The two-arc, two-node simple network: Performance of the fixed-point algorithm for 
BR-DUE, when the termination threshold e = 0.001. 
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e = 


e = 5 


£ = 10 


e = 15 


e = 20 


£ = 25 


Iteration 


18 


18 


16 


15 


14 


13 


Time (s) 


9.99 


9.63 


9.02 


8.45 


7.98 


7.51 



Table 2: The six-arc, five-node small network: Performance of the fixed-point algorithm for 
BR-DUE, when the termination threshold e = 0.001. 





£ = 


£ = 5 


£ = 10 


£ = 15 


£ = 20 


£ = 25 


Iteration 


22 


13 


11 


13 


10 


12 


Time (s) 


12.50 


7.79 


6.76 


7.95 


6.35 


7.69 



Table 3: The 19-arc, 13-node medium network: Performance of the fixed-point algorithm for 
BR-DUE, when the termination threshold e = 0.001. 





£ = 


£ = 5 


£ = 10 


£=15 


£ = 20 


£ = 25 


Iteration 


24 


7 


5 


6 


5 


4 


Time (s) 


459.3 


154.1 


119.5 


147.6 


122.2 


104.5 



Table 4: The 76-arc, 24-node Sioux Falls network: Performance of the fixed-point algorithm 
for BR-DUE, when the termination threshold e = 0.01. 

8 Conclusion 

In this paper, we have presented an analytical framework, known as (differential) variational 
inequality, that is proven to be a powerful tool for analyzing and computing dynamic user 
equilibria. Two natural extensions of the classical notion of simultaneous route-and-departure 
choice dynamic user equilibrium (Friesz et al. 1993[ ) are discussed in this paper. The first 
one is DUE with elastic travel demands (E-DUE); the second one is DUE with bounded 
rationality (BR-DUE). Both models have been studied in the literature to some extent but 
are never analyzed using the mathematical paradigm of differential variational inequality. 

We present the variational inequality (VI) and differential variational inequality (DVI) 
formulations of both E-DUE and BR-DUE. In particular, the VI formulation for BR-DUE 
is original in the existing literature. Analysis regarding the existence of both E-DUE and 
BR-DUE is most facilitated by expressing them as a variational inequality in properly defined 
Hilbert spaces and by invoking the Brouwer's fixed point theorem. Efforts in this direction 



have been elaborated by two papers Han et al. (2013c) and Han et al. (2013d). The VI 



formulations allow the still emerging theory of differential variational inequality |Pang and 
Stewart (2008) to be applied for analyzing and computing user equilibrium solutions, with 



surprisingly simple and elegant analytical results. We have demonstrated in this paper that 
both E-DUE and BR-DUE can be computed in continuous time in the DVI framework. The 
fixed-point-problem formulations also allow us to conduct rigorous analysis on the convergence 
of algorithms. 

One important direction of future research is to establish convergence results for a more 
general class of delay operators. Existing analyses on the convergence issue, described by either 



this paper or by Friesz et al. (2011), rely on some sort of monotonicity of the effective delay 



operator. As shown by the numerical results, the BR-DUE enjoys an improved convergence 
with positive e. Thus studying the monotonicity of the new operator $ e may provide a way 
of rigorously analyzing convergence and shed light on the monotonicity of the effective path 
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delay operator \&. 
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